Numerical precision radiative corrections to the Dalitz plot of baryon 
semileptonic decays including the spin-momentum correlation of the 

decaying and emitted baryons 

J. J. Torres 
Escuela Superior de Computo del IPN, 
Apartado Postal 75-702, Mexico, D.F. 07738, Mexico 

Ruben Flores-Mendieta 
Instituto de Fisica, Universidad Autonoma de San Luis Potosi, 
Alvaro Obregon 64, Zona Centro, San Luis Potosi, S.L.P. 78000, Mexico 

M. Neri, A. Martinez 
Escuela Superior de Fisica y Matemdticas del IPN, 
Apartado Postal 75-702, Mexico, D.F. 07738, Mexico 

A. Garcia 

Departamento de Fisica, Centro de Investigacion y de Estudios Avanzados del IPN, 
Apartado Postal 14-740, Mexico, D.F. 07000, Mexico 
(Dated: February 2, 2008) 



1 



Abstract 

We calculate the radiative corrections to the angular correlation between the polarization of the 
decaying and the direction of the emitted spin one-half baryons in the semileptonic decay mode to order 
{a/Tr){q/Mi), where q is the momentum transfer and Mi is the mass of the decaying baryon. The final 
results are presented, first, with the triple integration of the bremsstrahlung photon ready to be performed 
numerically and, second, in an analytical form. A third presentation of our results in the form of numerical 
arrays of coefficients to be multiplied by the quadratic products of form factors is discussed. This latter 
may be the most practical one to use in Monte Carlo simulations. A series of crosschecks is performed. 
This paper is organized to make it accessible and reliable in the analysis of the Dahtz plot of precision 
experiments involving heavy quarks and is not compromised to fixing the form factors at predetermined 
values. It is assumed that the real photons are kinematically discriminated. Otherwise, our results have a 
general model-independent applicability. 

PACS numbers: 14.20.Lq, 13.30.Ce, 13.40.Ks 
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I. INTRODUCTION 



The radiative corrections (RC) to spin one-half baryon semileptonic decays face three levels of 
complications. Despite the important progress achieved in the understanding of the fundamental 
interactions with the Standard Model [1], no first principle calculation of these corrections is yet 
possible. RC become then committed to model dependence and, what is worse, experimental 
analyses which use these calculations become model dependent, too. The second level comes from 
the fact that RC depend on the process characteristics, such as charge assignment of the baryons, 
type of the emitted charged lepton, size of the momentum transfer q involved, and whether real 
photons can be experimentally discriminated or not. RC also depend on the observable which is 
to be measured. All this requires RC to be recalculated every time the process characteristics and 
the observables are changed. At the third level one finds complications of a practical nature. It 
turns out that the final results of RC calculations are rather very inefficient to use or are long and 
tedious to the point that their use becomes unreliable. Fortunately, all the above complications can 
be solved rather satisfactorily. 

Although the model dependence of the virtual RC cannot be eliminated, an extension to baryon 
semileptonic decays of an analysis of Sirlin [2] of these corrections in neutron beta decay shows 
that to orders («/7r)(g/Afi)° and (a/n)(q/Mi), where Mi is the mass of the decaying baryon, 
the corresponding model-dependence amounts to several constants. These constants can all be 
absorbed into the already present form factors of the weak decay vertex. In addition, the theorem 
of Low [3] in its presentation by Chew [4] can be used to show that to these two orders of 
approximation the bremsstrahlung RC depend only on the non-radiative form factors and on the 
static electromagnetic multipoles of the particles involved. Accordingly, no model-dependence is 
introduced in this other part of the RC. Within these orders of approximation it is then possible to 
obtain final expressions for the RC that can be used in model-independent experimental analyses. 
The price is that it will be the effective form factors (which may be indicated by putting a prime on 
them) that can be experimentally determined. The separation of the original form factors from the 
model dependence of RC is then a theoretical problem only. It is in this sense that the first level of 
complications is put under control. 

To deal with the second level one must make an effort to calculate RC in a way as general 
as possible and to be able to obtain results which can be used directly to obtain the final results 
of other possible baryon semileptonic decays. In a recent publication |P] we showed that of the 
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six allowed charge assignments to the baryons when heavy quarks are involved, namely, A 
BH-Vu B+rv, A+ BH+uu A^ B-l+uu A++ B+l+uu and A+ B++l-Vu 

it is necessary only to calculate the RC to the first two, to which we shall refer to as charged 
decaying baryon (CDB) and neutral decaying baryon (NDB), respectively. The RC to the other 
four cases are then obtained from these two. We also showed that this property is valid up to order 
{a/n)(q/Mi), when / = e^, /x^, and r"^ and for any observable of baryon semileptonic decays. 
The problem of calculating RC is reduced considerably this way, although it will still be necessary 
to recalculate for different observables and whether real photons are discriminated or not in the 
first two cases. 

The third level of complications has been dealt with by computing analytically the triple 
integrals over the real photon variables. A numerical calculation of these integrals makes the 
application of RC to a Monte Carlo simulation practically impossible, because every time the 
values of the kinematical variables are varied those integrals must be recalculated. The analytical 
form of RC solves this problem. However, the results are very long and tedious and the use of 
this latter form may become unreliable. To control this it is very important that the analytical 
result be well organized and that it be crosschecked with the triple numerical integration form. A 
successful crosscheck allows the user to gain confidence on the analytical result and on its feeding 
into a Monte Carlo simulation. It may still be convenient to find a third presentation of RC, which 
would make their use more practical. 

From the above discussion it is clear that the calculation of RC to baryon semileptonic decays 
must be done following a program (see Ref. |Q] and references therein). In previous publications 
„e obtained ,he RC to the Dalitz plot of unpolarized decaying baryons up to order (a/.)(,/M,) 
MlM- In Ref. |8] we calculated to order (a/7r)(g/Mi)° the RC to the Dalitz plot with the angular 
correlation Si • p2 when the initial baryon is polarized along Si and the final baryon in emitted 
along p2 

In the present paper we want to attain two goals. The first one is to continue with our program 
and to calculate to order (a/7i)(q/Mi) the RC to the differential angular correlation Si ■ p2. The 
second one is to present the RC in the form of numerical arrays which should be applied to the 
quadratic products of form factors that appear in the RC, up to order {a/Tx){q/ Mi). We shall cover 
both CDB and NDB cases. 

The ordering of the paper is as follows. In Sec.|ll|we present the results to order {a/Tx){q/Mi) 
for the virtual RC. In SecHnlwe give the results for the bremsstrahlung RC in the triple numerical 
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integration form and combine them with the virtual RC results to obtain our first main result. 
In addition, we give the corresponding fully analytical results. In Sec. |lVl we perform several 
crosschecks and compare with other published results. In Sec.|Vlwe proceed towards our second 
goal. The last Sec.|VI|is dedicated to a summary and to concluding remarks. 

In order not to obscure the physics we have moved to Appendices A, B, and C the very many 
algebraic expressions that appear in the analytical results. In this paper we exhibit only new 
expressions. However, previously published expressions are required in these results. We do not 
reproduce them here. Instead, we give all the necessary references so that the reader can identify 
them correctly. The text and these appendices are organized so as not to obscure the physics and 
to make accessible the use of our results. Performing the analytical integrals is long an tedious. 
In order to help the reader interested in checking our results we have introduced the Appendix 
D, where the previous and new integrals can be identified. In Sees. IV and V we provide several 
numerical tables with the purposes of illustration and, more importantly, of helping the user to 
check his numerical results with ours. 



II. VIRTUAL RADIATIVE CORRECTIONS 



Our first purpose in this section is to review our notation and conventions. Next we shall discuss 
the virtual RC to the Si ■ p2 angular correlation over the Dalitz plot to order {a/Tx){q/ Mi). The 
uncorrected transition amplitude Mq for the baryon semileptonic decays 
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Gv is the Fermi decay constant multiplied by the appropriate Cabibbo-Kobayashi-Maskawa factor 
111. A and B are spin one-half baryons, / is the charged lepton, and ui is the accompanying 
antineutrino or neutrino as the case may be. ua, ub, ui or vi, and Vy or u^, are their corresponding 
spinors. The weak interaction vertex is 
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where = 7^(1 + 75), = (l/2)(7^7^ -7^,7^), and 7^ and 75 are Dirac matrices, q = pi -p2 
is the four-momentum transfer, and fiiq"^) and gi{q'^) are the vector and axial- vector form factors, 
respectively. Each form factor is assumed to be real in this work. The four-momenta and masses 
of the particles involved in ([T]) are pi = {Ei, pi) and Mi, p2 = {E2, P2) and M2, / = (E, 1) and 
m, and p^, = (E^, p^) and m^, respectively. Our calculations will be specialized to the center- 
of-mass frame of A. In this frame, pi, p2, I, and pu will also represent the magnitudes of the 
corresponding three-momenta, no confusion will arise from this. The directions of these momenta 
will be indicated by a caret, e.g., p2. 

Our approach to virtual RC follows the procedure of Ref. It has been discussed extensively 
in our previous works (see Ref. Id]), so only a few salient facts will be repeated here. The virtual 
RC can be separated into a model-independent part M^, and into a model-dependent one which 
amounts to six constants. These latter can be absorbed into the corresponding form factors of 
(|3l), this is indicated by a prime on Mq. The RC in M^, are finite in the ultraviolet, contain the 
infrared cutoff, and are gauge invariant. We shall limit ourselves here to exhibit explicitly only the 
new contributions of order (a/n)(q/Mi) to the Dalitz plot with the Si ■ p2 correlation. However, 
previous results are needed in the complete result. We shall give ample referencing to them. 

The transition amplitude with virtual RC is 

Mv = M; + M„. (4) 

The calculation of all the integrals over the virtual photon four-momentum that appear in My have 
been performed already to order {a/7c){q/Mi) in Ref. [6] for the CDB case and in Ref. |Q] for the 
NDB case. The corresponding results are compactly expressed as 

M.. = — [Mo$. + Mj,,$:], (5) 

where i = C, N separates the CDB and NDB cases, respectively. The matrix element M^c and 
the explicit forms of and $^ are found in Eqs. (8), (6), and (7) of Ref. 1^, respectively. The 
corresponding ones of M^n, $7v, and $^ are found in Eqs. (6), (7), and (8) of Ref. once the 
identifications = 2Re and = 2Re m0' are made. 

The Dalitz plot with virtual radiative corrections is now obtained by leaving E and E2 as the 
relevant variables in the differential decay rate for process (HJ) and specializing the result to exhibit 
explicitly the angular correlation §1 ■ p2. After making the replacement ua{pi) — * ^{si)uAiPi) 
in My [where the spin projector is given by = (1 — 75/i)/2], squaring the resulting 



amplitude, and rearranging terms we obtain for the differential decay rate 

In I TT J J 



(6) 



In Eq. ® the first two terms within curly brackets correspond to the unpolarized Dalitz plot. For 
i = C they can be found in Ref. where A'q, B[, and B'^.^, correspond to Eqs. (10), (11), 
and (12), respectively, of this reference. For i = N the unpolarized Dalitz plot can be found in 
Ref. [7], where B'^^ corresponds to Eq. (15). The spin-dependent part of Eq. ^ was obtained to 
order {a/n)(q/Mi)^ in Ref. [8]. There in Eq. (19) one can find the full expression for Aq. To the 
next order of approximation, however, there appear the new contributions, namely. 



B'2 = Ep2Qe + ElyoQ^ 



(7) 



B'^^ = Ep2Qs + ElyoQg, (8) 

and 

B';,^ = Mip2Qm + MilyoQN9. (9) 

The phase space factor of Eq. ^ is dfl = {Gy /2)[dE2 dE dQ2dipi/ {27r)^]2Mi, the cosine 
yo of the angle between the directions of the emitted baryon and the charged lepton is yo = 
[{E^y — P — pl]/{2p2l), and the neutrino energy is, by energy conservation, E^ = Mi — E2 — E. 

The Qi in Eqs. ©-©are functions of quadratic products of the form factors. They are new and 
are listed in Appendix A. 



ni. BREMSSTRAHLUNG RADIATIVE CORRECTIONS AND FINAL RESULTS 

The radiative process that accompanies dU) is 

A^B^Vl^, (10) 

where the real photon 7 carries four- momentum /c = (a;, k) and the neutrino energy is now E^ = 
El - u. 

The Dalitz plot for this four-body decay covers the three-body region of ^ and extends over 
it by a region where both E,, and uj are always non-zero simultaneously. We shall refer to this 
extension as the four-body region. A detailed discussion of these two regions as well as explicit 
expressions of their boundaries in the {E, E2) plane are given in Ref. jsl. Even if experiments have 
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no provision to detect the real photons in JTOb . a precise measurement of E and E2 still allows to 
discriminate against photons belonging to the four-body region. We shall assume in this paper that 
this is the case and shall restrict our calculations to the three-body region of (fTDt . 

In order to establish our notation and conventions and to make the necessary connections with 
our previous work, we must briefly review the derivation of the bremsstrahlung differential decay 
rate. According to the Low theorem jsl, the amplitude for process (flOb with contributions of orders 
1/k and /c° depends only on the form factors of the non-radiative amplitude © and on the static 
electromagnetic multipoles of the particles involved. The model dependence included by the real 
photon appears in new form factors which vanish at least linearly with k. These latter contribute 
to orders (a/7r)(g/Mi)^ and higher to the differential decay rate. Thus, within the approximations 
of this paper this part of the RC is model independent. The transition amplitude consists of the 
sum of three terms, namely, 

MiB = MiB, + MiB, + MiB,. (11) 

As in Secini the subindex i = C, N is used to distinguish CDB and NDB cases, respectively. The 
detailed expressions of Mcbi, Mcsa' ^cb^ are found in Eqs. (18), (19), and (20), respectively, 
of Ref. fl and of M^^,, MjvBa, and Mnb, are found in Eqs. (21), (22), and (23) of Ref. fl, 
respectively. Using the projector in Eq. (fTTT) . squaring the matrix element, performing the 

trace calculations, inserting the appropriate phase space factor, and indicating the integrations over 
the photon variables, the differential decay rate can be compactly given as 

dTiB = dT'.B - dr\'l (12) 

The analytical result to order («/7r)(g/Afi) of the unpolarized decay rate dV'^B was calculated in 
Ref. [6] for the CDB case and in Ref. ||7|] for the NDB case. They can be found in Eqs. (48) and 
(54) of such references, respectively. The polarized decay rate dV^^B was calculated analytically to 
order («/7r)(g/Mi)° in Ref. [SJ. The final results are given in Eq. (101) of this reference, for both 
CDB and NDB cases. 

The calculation to order {a/Tx){q/ Mi) of dV^^ is new. Let us now proceed with it. This decay 
rate consists of the sum of three terms 

dV^^ = dl^'iSi + dr^iu + drl%,,. (13) 

dV^^^j comes from the product Mj^^^Mj^i and it contains the infrared divergence and_the finite 
terms that accompany it. To extract them we follow the procedure used in Ref. L^, which 



extended the formalism introduced in Ref. fy] for Ki^ decays. The second and third summands in 
Eq. (fT3t come from of the product [M -^j^ + M -^^glMiSj for j = 2, 3 respectively. They are infrared 

(s) — 

convergent and are computed with standard techniques. The product M^-^gMj^s is left out because 
it contributes to orders («/7r)(g/Mi)^ and higher. The upper index s indicates where the 75 /i of 
S(si) is contained in these amplitudes. 

An important remark is in order here. It turns out that trying to compute the terms of 
order (a;/7r)(g/Mi) only and then adding them to the results of Ref. [s] is long and more 
cumbersome than doing from the start the full calculation containing both (a/n)(q/Mi)^ and 
{a/7i){q/Mi) contributions. Accordingly, our new expressions will contain the previous and the 
new contributions. It is then easy to verify that by eliminating the (a/7r)(g/Mi) terms in the new 
expressions one obtains the ones of Ref. [8]. 

The procedure to calculate the CDB and NDB cases differ substantially. We shall deal with 
them successively in the next two subsections 



A. Charged decaying baryon case 



The polarized radiative differential decay rate can be cast into the form 
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dnsi-f>2{B'2lco + C 



(14) 



Icq contains the infrared divergence and the finite terms that accompany it. It was calculated 
already, its explicit form is found in Eq. (52) of Ref. [8]. B'2 contains new (q/Mi) contributions. 
It coincides with Eq. (7) of the virtual RC. C^/^ consists of the sum of three terms, namely. 
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The integrations over the photon three-momentum are to be performed through the variables y = 
P2 • i, x = i • k, and the azimuthal angle Lpk of k. The traces of the square of the matrix elements 
give 
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[-k ■ p2(p2 + 2ly + 2cjk ■ P2) + lx]R2 + 
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k ■ P2(P2 + ly + ujk- P2)-R2 + 



UJ 



2Ep2 



Dp2 + 2ly{D - k • P2) 



Q 
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Dp2Uj E^^ . 

2Mi " 2MiE^^^ ~ 2 ■ P^^^^^^ 
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iP2 + ly + Ujk- P2)Ql8 



{E^ + (31 + pp2y + (3ujx)Q 
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Ml 
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k - P2 

I- I3x 



{(3E^ - P2y - I - ujx) + y 
{(3E^ +P2y + l + ujx) + y 
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/ (3y[x{El-D)+p2y + 
Ml 1 - /3x 

u; k ■ P2 

2Mi 1 - Px 
u 1 



21 



{E^-D + pl + I3p2y + pujx)Q22 
k ■ p2{Eu -pi- (3p2y - (3ux) + (3y{D - 2E,) 



2Mi 1 - /?x L 

+ ^TT^k ■ P2Q24 - 7TT^(P2 + + CUk ■ P2)Q25- 



23 



(19) 



2Mi ' 2Mi 

Here p = l/E, = -1 + - x^)/{l - /3x) + u/E, R2 = -l + {l- P'^)/{1 - (3x) - u/E, 
D = El + {\ + p2)- k, and uj = E/{2D), with F = 2^2/(1/0 - v)- 

The fomi factors of the vertex Q are contained in the Qi coefficients. These are collected in 
Appendix IaI 

The complete differential decay rate, containing the Dalitz plot with virtual and bremsstrahlung 
RC to order {a/n){q/Mi), is compactly expressed as 



dTc = dTcv + dT 



CB, 



(20) 



where the detailed expressions of dTcv and dVcs, containing si, can be traced starting at Eqs. Q 
and (fT4b . One can check that the infrared cutoff A contained in the virtual RC is canceled by 
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its counterpart in the bremsstrahlung RC. Eq. (l20b is model-independent to order («/7r)(g/Mi). 
The photon triple integrals of Eq. ([T6b remain to be performed numerically. This is our first main 
result, in the sense that it can already be used in a Monte Carlo simulation. It complies with 
all the requirements discussed in the Introduction to solve the difficulties of the first two levels. 
However, it still presents problems of the third level. The triple numerical integration form is 
still unpractical. This difficulty can be substantially solved because such triple integrations can be 
calculated analytically. 

We shall now proceed to obtain the analytical counterpart of the Si • p2 correlation contained 
in Eq. (l20b . Within our approximations all the form factors are constant and can be factored out of 
the very many triple integrals. A convenient rearrangement of the Cr of Eq. (fT31) is 

8 

1=1 

15 

Cii = (22) 



1=6 

and 

25 

Cm = 5^gA+3. (23) 

i=16 

The Qi are the quadratic functions of the form factors listed in Appendix|Al The triple integrals are 
contained in the Aj. We shall not detail here their explicit form in terms of such integrals. We only 
give their final analytical expressions and collect them in Appendix |Bl Many of these integrals 
have been performed already in our previous work, although some are new. To help the reader 
interested in following our calculations in more detail, we have given in Appendix O the general 
form of the triple integrals and a guide to identify their analytical counterparts in our previous 
work. Only the results for the new integrals are explicitly given in this appendix. In organizing 
Eq. (1211) with one running index i it was necessary to introduce A2 = 0, because Qs does not 
appear in this equation. 

The completely analytical result for the Dalitz plot in the differential decay rate of CDB, 
Eq. (EUt . can be compactly written as 



dTc = dVt 

where 



01 

- A'o'si ■ P2 + - (6c77 - 6c7/Si ■ P2) 
TT 



(24) 



Qci = B[{^c + Ico) + 5^1$ c + C'a, (25) 
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and 

Gcii = 5^($c + Ico) + 5^2*c + Ca^- (26) 

In this last equation $c, $'(^, B2 and the same of Eq. Jco is the one of Eq. (fT4l) . 

and C^"*^ is given by the sum of Eqs. (l2Tt - (l23t . In Eq. (|24l) and are the ones of Eq. Q 
and the analytic form of Qci is found in Eq. (48) of Ref. [6]. One can check that when (o/Mi) 
contributions in &cii and 6c/ are neglected, one obtains the order (g/Mi)° result of Ref. [8j. In 
particular B'^, B'^.^, and C^'^ become A'^, A'^, and D3{pi + pg) + D4{p2 + pi) of Eqs. (20), (21), 
and (96), respectively, of this reference. Let us now proceed with the second case. 



B. Neutral decaying baryon case 

The calculation of dv'^^^ proceeds in two ways. One possibility is to perform a straight-forward 
calculation using the tools described in the previous section. Another possibility is to use the 
approach introduced in Ref. |Q] to deal with the convergent pieces of dT^B- AH the (a;/7r)(g/Mi) 
terms can be obtained using the approximation 



P2 ■ k Pi ■ k {pi ■ kY 
and this will allow us to incorporate all the terms of order {a/u:){q/ Mi) that arise from this 
ratio. The advantage of this second possibility is that all the convergent terms of the NDB 
case are then obtained from their counterparts for the CDB case up to a few additional terms. 
This approximation, however, cannot be used in the divergent terms and hence we need standard 
techniques to calculate them. In the first term of Eq. (fT3t the infrared divergence is handled as in 
Ref. ||7|] and afterwards the approximation (l27b is used. One gets 

dT%j = -dnsi- MB2lm + Cj + Ci'^). (28) 

71 

The next term in Eq. (fT3t can be arranged using (l27t into 

dr^Nmi = -dnsi- MCii + Cjf). (29) 

71 

To calculate the third term in (fTSl) we can use the approximations W\ ^ 7a (/i +^7275) and p2 — pi. 
The traces that arise are practically the same as in dV^Q^jjj. However, there is a difference to order 
{a/7r){q/Mi) which gives rise to a Cjfj summand. Thus, one gets 

d^NBiii = -dnsi-p2 {Cm + C-S) . (30) 

71 
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The polarized decay rate becomes 



dV 



NB 



71 



NAJ- 



(31) 



Im contains the infrared divergence and the finite terms that accompany it. It was calculated 
already, its explicit form is found in Eq. (40) of Ref. BTJ. As before in Eq. ©, i?2 contains the 
/7r)(g/Mi) contributions. C^^^ is the same of Eq. (fT4ll of the CDB case and is defined as 



(a 



(32) 



After some tedious but straight-forward trace calculations their explicit forms are obtained. 
Ci'^^ becomes 



with i = I, II, III. Here D3 = 2{-gf + f\gi) and = 2{gf + fiQi). pi and p- are 
P2I 
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(38) 



(39) 
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The complete differential decay rate that contains the Dalitz plot of the NDB case including the 
Si ■ p2 correlation can be expressed compactly as 



dTN = dT]\fv + dTNB, (40) 

where the detailed expressions of dTjyv and dT^B containing si are traced starting at Eqs. Q 
and (BTT ). This is our second main result and the discussion of the previous subsection applies to 
it: its triple numerical integration form is still unpractical. This difficulty is solved by performing 
analytically the triple integrals contained in Eq. dSTT) . We have to concentrate only on C]^\ of this 
equation, all other terms in it have been given analytically already. Then, Eq. (l32b can be cast into 
the compact form 

dj^'^^ = D^pm + ApiV4, (41) 

where 

Pm = Pi + Pi I + Piii, (42) 

and 

Pn4 = p'i + p'ii + p'ni- (43) 

The explicit analytical expressions of the pi and p- are collected in Appendix ICl 

The completely analytical result for the Dalitz plot in Eq. (l40b can be put in parallel with 
Eq. (l24ll . namely. 



dr^v = dvt 



d 

- AqSi ■ P2 + - {QnI - ©AfZ/Sl ■ P2) 
TT 



(44) 



where 



and 



e^, = B[{^N + Im) + B'^^^'^ + C'^ + C^vA, (45) 



Qnii = 5^(<f TV + Im) + B';,^<^'^ + + (46) 

In this last equation ^tv, ^'n^ ^2 ^'^d ^7V2 the same of Eq. Q, C^^^ and J7V0 are the ones 
of Eq. (EB, and C^^\ is given in Eq. (1411) . In Eq. (1441) and are the ones of Eq. ^ and the 
analytical form of Qni is found in Eq. (54) of Ref. [7]. One can check that when {a / u:){q / Mi) 
contributions in Eq. (l44l) are neglected one obtains the (a/7r)(g/Mi)° result of Ref. [8]. In 
particular becomes A2 of this reference and C^]^ becomes zero. 
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IV. CROSSCHECKS 



There are several points we want to make in this section. One is that the analytical results 
are so long that it is important to check them. Another one is that there are some results already 
available in the literature |3] and we should compare with them. An even more important point 
is to provide the reader interested in using our result with numbers to be reproduced. 

To crosscheck the analytical results we use the triple numerical integration form of the RC. 
We make numerical comparisons of both forms by fixing the values of several form factors and 
of the Dalitz plot variables E and E2. A complete crosscheck requires the use of several choices 
of non-zero values for all the six form factors and a range of values of the pair (E, E2) over the 
Dalitz plot. Also, the comparison with the numerical results of Ref. Iiol should be made in the 
several cases covered there. All these crosschecks and comparisons were satisfactory and it is not 
necessary to display all the details here. Accordingly, we shall present a minimum of numerical 
tables and limit our discussion to them. 

For definiteness, we shall work with the decays E~ neV and A peV as examples of CDB 
and NDB cases. The reason for this is that numerical RC for these two decays were produced in 
Ref. 1 10]. We shall accordingly fix the form factors at the values used in this reference, namely, 
gi/fi = -0.34, /2//1 = -0.97 for ^ neV decay and gJf^ = 0.72, /2//1 = 0.97 for 
A peV. In addition, we use /i = 1 in ^ neV and /i = 1.2366 in A ^ peV and to compare 
with Ref. Iioi in both these decays we put g2 = and neglect and contributions. The values 
of the masses come from Ref. yj]. The anomalous magnetic moments of the baryons appear in 
our expressions of the RC. We use k{E-) = 0.376AMn, k{A) = 0.6130MAr, K{n) = 1.9130M^, 
and k{p) = — 1.7928MAr, where is the nuclear magneton. These values are extracted from the 
corresponding total magnetic moments reported in U] using Eq. (22) of Ref. t^. We neglected 
the anomalous magnetic moment of the electron, due to its smallness. 

As an example of the numerical crosscheck we display Table I for — > neV, where for 
generality we allowed g2,g3, fs 7^ 0. In the upper entries (a) we use the triple numerical integration 

(s) 

form to obtain the RC for C\ of the Si ■ p2 correlation covering a lattice of points over the Dalitz 
plot. The energies E and E2 enter through 6 = E/E^ and a = E2/M1. Em, a™"^, and cr™" are 
determined using the boundaries of the three body region given in Ref. jsl. The lower entries (b) 

is] 

contain the RC for the same C\ calculated with the analytical form. An inspection of this table 
shows an agreement to two decimal places and the third one being close. 
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(a) 



0.8077 


-0.0744 


-0.0811 


-0.0587 


-0.0256 


0.0101 


0.0420 


0.0650 


0.0741 0.0648 


0.0309 


0.8056 


-0.1350 


-0.1435 


-0.1048 


-0.0507 


0.0049 


0.0528 


0.0852 


0.0955 0.0778 


0.0253 


0.8035 




-0.1673 


-0.1294 


-0.0708 


-0.0090 


0.0443 


0.0802 


0.0910 0.0700 0.0102 


0.8014 




-0.1875 


-0.1534 


-0.0921 


-0.0254 


0.0327 


0.0718 


0.0831 0.0590 




0.7993 




-0.2065 


-0.1784 


-0.1152 


-0.0438 


0.0193 


0.0617 


0.0734 0.0459 




0.7972 






-0.2054 


-0.1407 


-0.0643 


0.0042 


0.0503 


0.0623 0.0306 




0.7951 






-0.2358 


-0.1697 


-0.0873 


-0.0124 


0.0379 


0.0499 0.0128 




0.7930 






-0.2720 


-0.2042 


-0.1141 


-0.0309 


0.0246 


0.0363 




0.7909 








-0.2480 


-0.1469 


-0.0520 


0.0106 


0.0211 




0.7888 








-0.3115 


-0.1913 


-0.0772 


-0.0035 


0.0041 




0.7867 










-0.2684 


-0.1110 


-0.0135 
















(b) 










0.8077 


-0.0744 


-0.0810 


-0.0587 


-0.0256 


0.0098 


0.0412 


0.0639 


0.0732 0.0643 


0.0309 


0.8056 


-0.1350 


-0.1435 


-0.1048 


-0.0507 


0.0047 


0.0521 


0.0842 


0.0947 0.0774 


0.0253 


0.8035 




-0.1673 


-0.1293 


-0.0706 


-0.0092 


0.0437 


0.0793 


0.0902 0.0697 


0.0102 


0.8014 




-0.1875 


-0.1533 


-0.0919 


-0.0255 


0.0322 


0.0710 0.0823 0.0587 




0.7993 




-0.2065 


-0.1783 


-0.1149 


-0.0437 


0.0189 


0.0609 


0.0727 0.0457 




0.7972 






-0.2053 


-0.1404 


-0.0640 


0.0040 


0.0497 


0.0617 0.0305 




0.7951 






-0.2357 


-0.1693 


-0.0869 


-0.0124 


0.0374 


0.0495 0.0127 




0.7930 






-0.2719 


-0.2036 


-0.1135 


-0.0308 


0.0243 


0.0360 




0.7909 








-0.2473 


-0.1460 


-0.0516 


0.0104 0.0210 




0.7888 








-0.3106 


-0.1900 


-0.0765 


-0.0034 0.0041 




0.7867 










-0.2666 


-0.1100 


-0.0134 






<5 


0.0500 


0.1500 


0.2500 


0.3500 


0.4500 


0.5500 


0.6500 0.7500 0.8500 0.9500 


^max 


0.8078 


0.8078 


0.8078 


0.8078 


0.8078 


0.8078 


0.8078 


0.8078 0.8078 


0.8078 




0.8043 


0.7978 


0.7925 


0.7884 


0.7857 


0.7847 


0.7854 0.7884 0.7939 


0.8023 



TABLE I: Values of in S 



-> nev decay by (a) integrating it numerically and (b) evaluating it 
analytically. C^f^ is given in units of GeV^. The form factors have been given the arbitrary values /i = 1.0, 
/2 = -0.97, h = -0.778, gi = -0.34, 32 = 0.987, and 53 = -1.563. 



To proceed with the comparison with Ref. [101 we must use the difference defined there, 
namely, 

5aB{E, E2) = aB{E, E2) - ao{E, E2), (47) 



where 



aBiiE,E2) 



A'^{E,E2) + {a/7i)Qui{E,E2] 



(48) 



A',{E,E2) + ia/7T)QuiE,E2) 
and as before i = C,N and ao{E, E2) = -Aq{E, E2) /A'q{E, E2) for both values of i. 

One may interpret aBi{E, E2) as the asymmetry parameter of the emitted baryon at (£", E2 
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points of the Dalitz plot. Here we must choose the same {E, E2) points as in Ref. [^1 and, as 
already mentioned, use the same values of the form factors. However, it should be stressed that 
our final results are not compromised to fixing the values of the form factors. 

Before proceeding with a detailed comparison with the numbers of this reference, there is 
a point that must be kept in mind. The approximations used in our work and in Ref. 1 10] are 
not quite the same. We used the Low theorem to calculate the bremsstrahlung RC and in this 
reference it was assumed that both baryons involved were point-like and higher (a;/7r)(g/Mi)"' 
contributions (n > 2) were included in this part. Another interesting thing is to compare our order 
(a/7r)(g/Mi) results with our previous order {a/n){q/ Mif results. As explained earlier these 
latter are reproduced here when the (a/7r)(g/Mi) contributions are neglected. 

We performed many comparisons and, as before, there is no need to present all the details. 
One example, the A peV case is enough for this discussion. The results are displayed in Table 
II. In the upper (a) part only the order (a;/7r)(g/Afi)° is given. Both this order and the order 
(a/7r)(g/Mi) contributions are added in the middle part (b). The numerical results of Ref. \\\ 
are reproduced in the lower part (c). A numerical crosscheck was also performed in producing 
parts (a) and (b). We do not reproduce it here, the agreement was as good as in Table I. 

An inspection of Table II shows that the order {a/Tx){q/ Mi) is systematically perceptible at the 
second significant digit and even at the first one. In comparing with Ref. flO^, one can see a better 
agreement with the middle table (b). The agreement at the first significant digit improves as the 
RC grow in size and also the variations in the second digit become smaller. There are differences, 
however. They may be explained as due to the different approximations used. Also, comparing 
entries (a) and (b) one may conclude that for light quark hyperon semileptonic decays the order 
(a/7r)(g/Mi) is perceptible enough and that when heavy quarks are involved contributions of this 
order become relevant in precision experiments. 

Let us now turn to a different form to use our results. A form which may provide a more 
efficient use of them in a Monte Carlo simulation and which still is not compromised to fixing 
values of the form factors, as was the case in Table n. 

V. NUMERICAL FORM OF THE RADLVTIVE CORRECTIONS 

We now come to our second goal in this paper. In the previous sections we have obtained the 
RC to CDB and NDB in two forms. The first one has triple integrals over the real photon variables 
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a (a) 



0.8530 


0.1 


0.0 


0.0 


0.0 


0.0 


0.1 


0.1 


0.2 


0.4 


0.8518 


1.2 


0.2 


0.1 


0.1 


0.2 


0.2 


0.3 


0.5 


1.0 


0.8505 




0.5 


0.3 


0.3 


0.3 


0.4 


0.6 


0.9 


1.5 


0.8492 




1.0 


0.5 


0.5 


0.5 


0.6 


0.8 


1.2 


1.6 


0.8480 




1.7 


0.9 


0.7 


0.7 


0.8 


1.1 


1.4 


1.4 


0.8467 






1.3 


1.0 


1.0 


1.1 


1.3 


1.5 


0.5 


0.8454 






1.9 


1.4 


1.3 


1.4 


1.5 


1.5 




0.8442 








1.9 


1.7 


1.7 


1.7 


1.0 




0.8429 








2.7 


2.3 


2.1 


1.6 






0.8416 










3.4 


2.4 


0.1 
















(b) 










0.8530 


0.5 


0.2 


0.2 


0.2 


0.1 


0.1 


0.1 


0.2 


0.3 


0.8518 


2.0 


0.5 


0.4 


0.4 


0.4 


0.4 


0.5 


0.6 


0.9 


0.8505 




0.9 


0.7 


0.6 


0.6 


0.7 


0.8 


1.0 


1.3 


0.8492 




1.5 


1.0 


0.9 


0.9 


0.9 


1.1 


1.3 


1.3 


0.8480 




2.3 


1.3 


1.1 


1.1 


1.2 


1.3 


1.4 


1.0 


0.8467 






1.8 


1.5 


1.4 


1.5 


1.5 


1.5 


0.2 


0.8454 






2.3 


1.8 


1.7 


1.7 


1.7 


1.3 




0.8442 








2.2 


2.0 


1.9 


1.7 


0.7 




0.8429 








2.8 


2.4 


2.1 


1.5 






0.8416 










3.2 


2.2 


0.1 
















(c) 










0.8530 


0.2 


0.1 


0.1 


0.1 


0.1 


0.1 


0.1 


0.2 


0.3 


0.8518 


1.9 


0.4 


0.3 


0.3 


0.3 


0.4 


0.5 


0.6 


1.0 


0.8505 




0.8 


0.6 


0.5 


0.6 


0.6 


0.8 


1.0 


1.3 


8493 




1 4 


Q 


8 


u.o 


Q 


\ \ 


1 % 


1 4 


0.8480 




2.4 


1.3 


1.1 


1.1 


1.2 


1.3 


1.5 


1.1 


0.8467 






1.7 


1.4 


1.4 


1.5 


1.6 


1.5 


0.3 


0.8455 






2.3 


1.8 


1.7 


1.7 


1.7 


1.3 




0.8442 








2.2 


2.0 


2.0 


1.7 


0.8 




0.8429 








2.8 


2.5 


2.2 


1.5 






0.8417 










3.2 


2.3 


0.1 







<5 0.0500 0.1500 0.2500 0.3500 0.4500 0.5500 0.6500 0.7500 0.8500 0.9500 

^max 0.8536 0.8536 0.8536 0.8536 0.8536 0.8536 0.8536 0.8536 0.8536 0.8536 
^mm 0.8516 0.8479 0.8450 0.8428 0.8414 0.8410 0.8416 0.8433 0.8464 0.8508 

TABLE II: 100(5qb(£', E2) with RC over the three-body region in A — > peV decay, (a) gives the RC to 
order ia/7r)(g/Mi)0, (b) gives the RC to order (a/7r)(g/Mi), and (c) corresponds to the RC computed in 
Ref. Ildl. 
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ready to be performed numerically. The second one is fully analytical. Although this latter one is 
already practical it is still long and tedious. It still requires that the RC be calculated within the 
Monte Carlo simulation every time E and E2 are varied. This is much faster than performing the 
triple integrals, but, it still represents a non-negligible computer effort. We shall now discuss a 
third form of the RC that may be more practical to use. 

For fixed values of E and E2, Eqs. (El) and for the CDB case and Eqs. (|45l) and for 
the NDB case take the form 

6 

= E (49) 

i<j=i 

because they are quadratic in the form factors. For the sake of simplify, in Eq. ( l49l) we have 
momentarily redefined gi = fi,g2 = /s, and = /q. Notice that the restriction i < j reduces 
the sum in Eq. (l49l) to 21 terms. The index m takes the values m = CI, CII, NI, and NIL The 
third form of RC we propose consists of calculating arrays of the a™ coefficients determined at 
fixed values of (E, E2) and that these pairs of {E, E2) cover a lattice of points on the Dalitz plot. 

To calculate the coefficients it is not necessary to rearrange our final results, either analytical 
or to be integrated, so that they take the form ( l49l) . One can calculate them following a systematic 
procedure. One chooses fixed {E, E2) points. Then one fixes /i = 1 and /j = 0, i 7^ 1 and obtains 
a^, one repeats this calculation for /2 = 1, = 0, z 7^ 2 to obtain a^, and again until /e = 1, 
/j = 0, i 7^ 6 and is obtained. Next, one repeats the calculation with /i = 1, /2 = 1, fi = 0, 
i ^ 1,2 and from this results one subtracts and a^, this way one obtains the coefficient a^. 
One repeats this last step changing i and j until all the interference coefficients a^, i 7^ j, have 
been calculated. 

To illustrate all this and to further discuss it we have produced arrays presented in two tables, 
selecting in each one ten points {E, E2) over the Dalitz plot. We have chosen two examples, 
A peV of a NDB case which is displayed in Table III and A+ Ae+z/ of a CDB case which 
is displayed in Table IV. This latter also serves as an example of a heavy quark decay. As in the 
previous section, the more important purpose is to provide the user with numbers to compare with. 
The arrays of these two tables were obtained using the RC in the analytical form. In the A+ case 
we used the formulas for the charge assignments A^, B^, of the CDB case of the previous 
sections and then applied the rules of Ref. jsl to obtain the results for the charge assignment , 
B^, /+ of this particular case. In these tables we have restored our standard notation for the axial- 
vector form factors gi, g2, and g^. The masses used are those of Sec. IIVI Mi(A+) comes from 
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Ref. fl], and we assume an estimate for /t(A+) = O.IIOGMat. 

The first fact tliat appears in these tables is that the RC do not depend on the form factor 
products /|, gl, /1/3, /s/s, QiQz, 9293, fi93, f293, and /g^g in Tables III and IV. The non- 
appearance of these products cannot be seen easily in our final results of Sec. ITTD The other 
fact is that the non-zero RC to each form factor product vary appreciably from one {E, E2) point 
to another. This means that replacing the precision results of Sec. [nil with an array of only a few 
columns over the Dalitz plot is far from satisfactory. Therefore the lattice of (E, E2) points must 
be much finer than only a few points. 

The use of this third presentation of RC is very practical in the sense that such RC can be 
calculated separately and only the arrays should be fed into the Monte Carlo simulation. However, 
in a precision experiment possibly involving 150, 200, and even 300 bins over the Dalitz plot 
the number of columns in the RC arrays should be at least just as many. It may be required that 
several columns be produced in finer subdivisions within each bin, possibly 4, 8, or even more. For 
example, one may require that the numerical changes of the coefficients between neighboring 
[E, E2) points do not exceed two decimal places within rounding of the third decimal place. 

To close this section let us stress that none of the three forms of our RC results is compromised 
to fixing from the outset values for the form factors when RC are applied in a Monte Carlo 
simulation. To fix them at prescribed values may be not too bad an assumption for hyperon 
semileptonic decays, but it is not acceptable at all for decays involving heavy quarks where the 
Cabibbo theory yj] is no longer reliable for fixing the form factors. 

VI. SUMMARY AND CONCLUSIONS 

We have obtained in Sees. |ll| and |lll| the RC to the angular correlation §1 ■ p2 to order 
{a/'K){q/ Ml). Our final results are given in two forms. The first one is the triple numerical 
integration form, in which the integrations over the real photon variables are explicitly exhibited 
and remain to be performed numerically. The second one is the analytical form where those 
integrations have all been calculated analytically. We covered two cases, the CDB and the NDB 
ones whose final results are given in Eqs. (l24b and (l44b . respectively. The analytical results are 
very long and tedious. To make their use more accessible we have collected the numerous Qi and 
Aj algebraic expressions which appear in the CDB case in Appendices IaI and IbI respectively. The 
NDB case uses these expressions and also the pi ones. These latter were collected in Appendix O 
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(0.05,0.8518) 


(0.35,0.8518) 


(0.65,0.8518) 


(0.95,0.8518) 


(0.25,0.8480) 


(0.55,0.8480) 


(0.75,0.8480) 


(0.45,0.8442) 


(0.65,0.8442) 


(0.55,0.8416) 


f! 


6.812 


X 10-4 


1.550 


X 10-4 


-1.024 


X 10-4 


6.505 


X 10-4 


7.218 


X 10-4 


-1.622 


X 10-4 


-6.987 


X 10-5 


7.198 


X 10-5 


-1.227 


X 10-4 


-4.887 


X 10-5 


fl 


1.604 


X 10"^ 


-1.875 


X 10-4 


-5.732 


X 10-4 


2.204 


X 10-3 


6.787 


X 10-4 


-1.268 


X 10-3 


-3.588 


X 10-4 


-7.888 


X 10-4 


-5.807 


X 10-4 


-2.633 


X 10-4 


fi 




0.000 




0.000 




0.000 




0.000 




0.000 




0.000 




0.000 




0.000 




0.000 




0.000 


9l 


7.470 


X 10-2 


-8.900 


X 10-3 


-2.635 


X 10-2 


1.026 


X 10-1 


4.812 


X 10-2 


-9.116 


X 10-2 


-2.537 


X 10-2 


-1.250 


X 10-1 


-9.181 


X 10-2 


-2.073 


X 10-1 


9l 


1.882 


X 10-=^ 


-2.199 


X 10-4 


-6.722 


X 10-4 


2.585 


X 10-3 


1.227 


X 10-3 


-2.293 


X 10-3 


-6.488 


X 10-4 


-3.131 


X 10-3 


-2.305 


X 10-3 


-5.219 


X 10-3 


9l 




0.000 




0.000 




0.000 




0.000 




0.000 




0.000 




0.000 




0.000 




0.000 




0.000 


hh 


2.010 


X 10-^ 


-4.892 


X 10-5 


-5.756 


X 10-4 


2.395 


X 10-3 


1.258 


X 10-3 


-1.173 


X 10-3 


-3.495 


X 10-4 


-4.985 


X 10-4 


-5.694 


X 10-4 


-2.477 


X 10-4 


/1/3 




0.000 




0.000 




0.000 




0.000 




0.000 




0.000 




0.000 




0.000 




0.000 




0.000 


f2h 




0.000 




0.000 




0.000 




0.000 




0.000 




0.000 




0.000 




0.000 




0.000 




0.000 


9192 


-2.362 


X 10-2 


2.889 


X 10-=* 


8.504 


X 10-3 


-3.252 


X 10-2 


-1.509 


X 10-2 


2.906 


X 10-2 


8.219 


X 10-3 


3.975 


X 10-2 


2.911 


X 10-2 


6.574 


X 10-2 


9i9-i 




0.000 




0.000 




0.000 




0.000 




0.000 




0.000 




0.000 




0.000 




0.000 




0.000 


9293 




0.000 




0.000 




0.000 




0.000 




0.000 




0.000 




0.000 




0.000 




0.000 




0.000 


fi9i 


1.683 


X 10-1 


3.604 


X 10-1 


1.875 


X 10-2 


-5.329 


X 10-2 


2.878 


X 10-1 


8.435 


X 10-2 


-1.250 


X 10-1 


1.258 


X 10-1 


-8.629 


X 10-2 


-3.368 


X 10-2 


/192 


-1.763 


X 10-=* 


-6.857 


X 10-=* 


6.870 


X 10-4 


-9.573 


X 10-4 


-1.808 


X 10-2 


-4.425 


X 10-3 


9.176 


X 10-3 


-1.410 


X 10-2 


1.046 


X 10-2 


5.085 


X 10-3 


fi9-i 




0.000 




0.000 




0.000 




0.000 




0.000 




0.000 




0.000 




0.000 




0.000 




0.000 


/291 


-2.482 


X 10-2 


-8.236 


X 10-3 


-6.172 


X 10-3 


2.989 


X 10-2 


-2.846 


X 10-2 


-1.758 


X 10-2 


-4.052 


X 10-3 


-1.703 


X 10-2 


-9.910 


X 10-3 


-5.365 


X 10-3 


/292 


3.970 


X 10-^ 


1.329 


X 10-3 


9.192 


X 10-4 


-4.778 


X 10-3 


4.634 


X 10-3 


2.719 


X 10-3 


4.868 


X 10-4 


2.659 


X 10-3 


1.401 


X 10-3 


7.390 


X 10-4 


/293 




0.000 




0.000 




0.000 




0.000 




0.000 




0.000 




0.000 




0.000 




0.000 




0.000 


h9i 


-5.014 


X 10-5 


-9.348 


X 10-5 


-6.324 


X 10-5 


-5.962 


X 10-^ 


-1.323 


X 10-4 


-1.003 


X 10-4 


-2.914 


X 10-5 


-9.289 


X 10-5 


-2.412 


X 10-5 


-9.643 


X lO-'' 


h92 


-4.271 


X 10-5 


-8.147 


X 10-5 


-5.522 


X 10-5 


-5.513 


X 10-^ 


-1.136 


X 10-4 


-8.644 


X 10-5 


-2.517 


X 10-5 


-7.876 


X 10-5 


-2.061 


X 10-5 


-8.140 


X 10-'^ 


f393 




0.000 




0.000 




0.000 




0.000 




0.000 




0.000 




0.000 




0.000 




0.000 




0.000 



TABLE III: Numerical arrays of the coefficients afj^^ in GeV^ of Eq. d49l evaluated at ten points {E, E2) (headings of columns) over the polarized 
Dalitz plot of A — > pe'v decay. 



(0.15,0.5995) (0.45,0.5995) (0.75,0.5995) (0.95,0.5995) (0.25,0.5602) (0.55,0.5602) (0.85,0.5602) (0.45,0.5210) (0.75,0.5210) (0.65,0.4948) 



n 


-1.532x10" 


-2 


1.036x10" 


'1 


6.859x10" 


-3 


-3.894x10" 


-1 


7.582x10" 


-2 


2.298x10" 


1 


-3.036x10" 


-1 


2.295x10" 


- 1 


-3.851x10" 


-2 


6.446x10-2 


/I 


-4.818x10" 


-3 


2.436x10" 


-1 


1.849x10" 


-2 


-8.624x10" 


1 


2.289x10" 


-1 


5.389x10" 


-1 


-6.710x10" 


-1 


5.611x10" 


- 1 


-8.135x10" 


-2 


1.481x10"^ 


/I 


-0.000 




-0.000 




-0.000 




-0.000 




-0.000 




-0.000 




-0.000 




-0.000 




-0.000 




-0.000 


9! 


-1.466x10" 


-2 


1.076 




8.152x10" 


-2 


-3.808 




1.534 




3.546 




-4.416 




7.830 




-1.139 




1.004x10+1 


9l 


-5.571x10" 


-3 


2.817x10" 


'1 


2.138x10" 


-2 


-9.972x10" 


-1 


3.943x10" 


-1 


9.282x10" 


1 


-1.156 




2.047 




-2.967x10" 


-1 


2.632 


9l 


-0.000 




-0.000 




-0.000 




-0.000 




-0.000 




-0.000 




-0.000 




-0.000 




-0.000 




-0.000 


hh 


-3.445x10" 


-2 


3.133x10" 


-1 


2.160x10" 


-2 


-1.159 




2.497x10" 


-1 


6.956x10" 


-1 


-9.031x10" 


-1 


7.047x10" 


- 1 


-1.130x10" 


-1 


1.941x10"! 


hh 


0.000 




0.000 




0.000 




0.000 




0.000 




0.000 




0.000 




0.000 




0.000 




0.000 


hh 


0.000 




0.000 




0.000 




0.000 




0.000 




0.000 




0.000 




0.000 




0.000 




0.000 


9192 


8.296x10" 


-3 


-1.114 




-8.738x10" 


-2 


3.897 




-1.571 




-3.659 




4.514 




-8.037 




1.150 




-1.030x10+1 


9193 


0.000 




0.000 




0.000 




0.000 




0.000 




0.000 




0.000 




0.000 




0.000 




0.000 


9293 


0.000 




0.000 




0.000 




0.000 




0.000 




0.000 




0.000 




0.000 




0.000 




0.000 


fi9i 


9.235 




9.436x10" 


'1 


-9.728 




-2.410 




8.070 




-2.512 




-5.571 




5.378 




-6.121 




-2.391 


fl92 


-4.429x10" 


-1 


8.400x10" 


-2 


7.031x10" 


-1 


-8.759x10" 


-1 


-1.671 




7.110x10" 


■ 1 


5.766x10" 


-1 


-2.146 




2.184 




1.134 


fl93 


0.000 




0.000 




0.000 




0.000 




0.000 




0.000 




0.000 




0.000 




0.000 




0.000 


h9i 


-4.903x10" 


-1 


-2.339x10" 


-1 


1.770x10" 


-2 


3.505 




-2.921x10" 


'2 


-6.449x10" 


-1 


3.069 




1.706 




7.521x10" 


-1 


2.504x10-1 


h92 


3.501x10" 


-1 


1.071x10" 


-1 


-1.408x10" 


-1 


-1.853 




2.102x10" 


-1 


2.404x10" 


-1 


-1.777 




-7.944x10" 


-1 


-5.591x10" 


-1 


-1.778x10"! 


h93 


0.000 




0.000 




0.000 




0.000 




0.000 




0.000 




0.000 




0.000 




0.000 




0.000 


h9i 


-2.067x10" 


-2 


-2.617x10" 


-2 


-1.189x10" 


-2 


-4.642x10" 


-4 


-3.808x10" 


'2 


-3.174x10" 


-2 


-3.158x10" 


-3 


-3.269x10" 


-2 


-4.636x10" 


-3 


-3.500x10-3 


h92 


-1.261x10" 


-2 


-1.597x10" 


-2 


-7.253x10" 


-3 


-2.833x10" 


-4 


-2.147x10" 


-2 


-1.790x10" 


-2 


-1.781x10" 


-3 


-1.704x10" 


^2 


-2.417x10" 


-3 


-1.731x10-3 


h93 


0.000 




0.000 




0.000 




0.000 




0.000 




0.000 




0.000 




0.000 




0.000 




0.000 



TABLE IV: Numerical arrays of the coefficients afj^^ in GeV^ of Eq. ( l49t evaluated at ten points {E, E2) (headings of columns) over the polarized Dalitz 
plot of ke^v decay. 



Our analytical results were crosschecked and compared with other results available in literature. 
This we have done in Sec. HVj We have limited ourselves to discuss the decay neV as an 

example of the crosschecks and A peV as an example of the comparisons with Ref. liioi . 
In addition, in this latter decay we also included a comparison between our RC to orders 
(a/7r)(g/Mi)0 and (a/7r)(g/Mi). 

We have discussed in Sec. IVl another possibility to use our results in an experimental analysis. 
One can calculate the numerical factors of the quadratic products of the form factors that appear 
in the RC at fixed values of (E, E2). These factors can be organized in arrays to be multiplied 
upon such products, covering a lattice of {E, E2) points over the Dalitz plot. We discussed two 
examples of this possibility, a CDB one and NDB one. 

Apart from illustration purposes, the several tables in this paper provide numbers to compare 
with. Also, our calculations rely heavily on previous results. Apart from discussions in text, we 
have given in Appendix details to allow the identification of our previous and new analytical 
results for the many integrals. 

To close, let us recall that our results are general within our approximation. They can be applied 
in the other four charge assignments of baryons involving heavy quarks and whether the charged 
lepton is e^, js^, or t^. They are model independent and are not compromised to fixing the form 
factors at prescribed values. The above calculations should be extended to cover precision RC in 
the Si ■ i correlation and in the four body region [[l2l[l . We hope to return to these cases in the 
near future. 
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APPENDIX A: COLLECTION OF THE Q, COEFFICIENTS 



The coefficients Qi introduced in Sees. II and III are long quadratic functions of the form 
factors. The coefficients Qi, . . . ,Qj have been computed in previous works [sl and can be 
found there. The new coefficients are listed below. They read 



E2- M2- Pp2yo 



Ml 

— {F1F2 — G1G2 



E2 + M2- pp2yo 



f3p2yo 



+ F2G1 



Qi = Fl 



1 + 



M2 

Ml 



Ml 

M1 + M2 + EI-E 
Ml 

E2-M2 



— F1G2 



Ml 

Ml -M2 + EI-E 



E2 - pP2yo 



Ml 



E 



+ F1G2 



E2-M2 



Ml 

+ (F1F2 - G1G2) 



+ Gl 



E^. - E 



Qs = F^ 



Eo - Mo 



Ml 



— F2G2 



E'. 



Q9 = Fl 



Ml 
E2-M2 



+ Gl 
— F2G1 



E 

MiE 

E2 + M2 



2MiE_ 

\_M2 

Ml 
- F2G1 



Ml 
— F2G2 



2E^. 



2MiE 
2 



1 



Ml 2MiE 
E2 + M2 



E 

E2 + M2 
Ml 



- 2FiGi 

El-E - 
E 



El - E 
E 



+ 2F1G1 



Ml 



+ Gl 



Ml 

E- M1-M2 
Ml 

'E2 + M2' 



En 



Ml 
+ F3G1 



+ F1G2 



E - Mi + M2 



Ml 



E{E2 + M2) 



+ F2G1 



E2 + M2 



Q 



N8 



F'l 



Ml 

{Mi-E){E2-M2 



Ml 
F3G1 



El 

Ml 



+ 2F1G1 
1 



Eo 



Ml 
E2 + M2 



Ml 



— F1G2 



E2-M2 
Ml 



+ F1G2 
— F2G2 



Ml 

M2{E2- M2- El 



+ Gl 



Ml 

{Mi-E){E2 + M2] 



EoE?. 



Q 



N9 



-Fl 



Ml 

M2{E2 - M2 



Ml 

+ F3G1 



+ F2G1 



Ml 

E2{E2 



+ Gl 



Mo 



Ml 



+ F,Gi 



+ F2G1 

'E{E2 + M2) 
Ml 

M2{E2 + M2)' 

Ml 

' M2 + E2 
Ml 



Ml 

M2{E2 + M2-El) 
Ml 





Ml 


El- 


+ 2F1G1 






M 


Mi_ 



+ 2^16-1 

E2 



Ml 



— FIG2 



E2{E2 - M2) 
Ml 



1 - 



Ml 



24 



and 



Qio — —F2G1 + F1G2 + F2G2, 
'E2 + M2' 



Q 



11 



Ml 

Q12 — 2FiGi, 



G1F3, 



Q 



13 



-F^ 



E2 - M2 



-G\ 



— F1G2 



E 

E2-M2 
E 

Qu ■ 
Qib = '2F^ 



E2 + M2 



E 



+ 2FiGi 



Eo 



F3G1 



E2 
Ml 



E _ 

M2 + E2 
E 



+ F2G1 



E2 + M2 
E 



—Ff — G\ — F1F2 + G1G2, 



E2 - M2 
Ml 



+ 2Gl 



E2 + M2 
Ml 



<5l8 = 



Qie = fi{92 - gi) - f29i, 
Qn = fi92 + fsQi, 

l(fi-9!) + f2(fi + 9i)-gi(f3-g2), 



Ql9 = 2/igri 

Q20 = -2gl 



i Ki 

2Mi 

1 Ki 

2Mi 



Ml, 
Ml, 



Q21 = 7.(/i' - gl) + /2(/i - 51) + giih + g2), 



Q22 = ifi - 9i){f2 - 92) + Mi^(/i - - Mi^(/f - gf), 

e e 

Q23 = -(/i + gi)(f2 - 52) + Mi^(/i + - M,^{fl - gl), 

e e 



Q24 = -(/i - ^1)' + ^1(2/1 + 3/2 + 2/3 + ^2 - 2^1) - /i(/2 + g2) 
e 



Mr-{h - gi? + Mi^(/f - gf) - AM^gl 



Q2, = -if? - gl) - (/i + gi){f2 + g2) + 2^i(/3 - 72) 

+ Mi^(55i^ + fl + 2/1^1) - Mi^(/i^ - gl). 



Here, ki and K2 denote the anomalous magnetic moments of the decaying and emitted baryons, 
respectively. 
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The tildes on Qq and indicates that contributions of order [q/MiY and higher have been 
subtracted. Also, Qs, ■ ■ ■ , Q2b have only contributions up to order q/Mi. 

Although we have not made it explicit, in the above expressions the primed form factors, 

containing the model dependence of virtual RC should be used. This is valid to order (a/vr)^ 
rearrangements. In the coefficients Qq, ■ ■ ■ ,Qi5 we have used the Harrington's form factors 
Fi, Gi. They are related to the Dirac's form factors fi, gi as Fi — /i + (1 + M2/Mi)/2, 
Gi = ^1 - (1 - M2/Mi)g2, F2 = -2/2, G2 = -2^2, ^3^/2 + fs, ^3 = ^2 + ^3. 

APPENDIX B: COLLECTION OF THE FUNCTIONS 

Here we give the analytical expressions of the Aj functions that appear in Sec. HI in the 
analytical form of the RC to the polarized decay rate. 



Ai = -EW, 



A3 = ^ - 1)X12 + 2x11 - Xio] , 



Elpl 



.2 r 




A4 



2Mi 



As 



El 



pIy^ + 2Z2 + 



2p2l' 



2Mi 



E 




[p19o + {E + El)Zi - EZ2 - fp2Yi] , 



Ay 




A, 



■8 




A, 



9 




Aio — 2^Cii, 
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Aii = ^/3p2(7o-^^3), 



A 



12 



X 1 

-KlO - /?^3 - + ^{X21 - X20) + -^^2 



A 



13 



2Mi 



V ^2 ^21 -^ + -^° 



-^7 + E{e, - 62) 



Ai4 = 



2Mi 

P2^ 



^2 

7o - PW^ + — -r]o 



+ 



X,- 



+ 



Ml 



2 M^E 2Mi 
_|x, + |[/(yo-l)-2p2]], 



+ El 
AMiE 



pjOr + 2plE{94 - 9s) + 2(21 - 7X3 



+ 



[4(^2/0 +^2)^2 +P2X21] , 



I 



+ 



X4 P2I' 



2Mi 4Mi 



(y 



A _ ^ ^ 



A 



17 



M 

4Mi 



^ - 2£;77o + + E'^) [9^ + 2E{9, - 9s)]] , 



A18 = 4^ [^3 - 2£;°X2] , 



Ai9 - ^ [X20 - X21 + 2El{xn - Xio) + ^{(21 - 2S°Cii)] , 



A 



20 



E 



(3fp2(Y, - Y,) - {El + l(3){xii - Xio - /3Cii) 



+ (5{E + £;°)(Cii - Cio) + ^WP2{1 - 2/o)(^o + m) + /9'^P2/ 



/^/^(n + - + ^p^n + 1Z2 
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A 



22 



E 

Ml 



+ p{E - £;°)(Cii - Clo) + + fPP2{i - vl) 



A23 = — 
^ Ml 



^ + WE'. + l)Xii -iE + E'^){Cu - Cio) + lp2l{l - yl) 



-P2l(yo - 1)00 + ElCio - {El + Z/3)Cii - (3lpll 



A 



P2I 
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4Mi 

P2 



2z/oXii + ^(^0 + 2r7o)(l - t/o) + -(^ + ^° - /?P2yo)Cii 

P2 



- — Clo - + 131)^ + PiE + + TVi + 2/3/P2/ 



P2^ 



P2^ 



A26 = ^ [(^° - ^/?)X21 - P{El - E)U + P2/iVi - X31 + /3C3l] , 



and 



A27 = ^[^°X20+P2/iV2], 



A28 = -^A^2. 



4Mi 



In these Aj we introduced the definitions 

,2 



-^2 = "^^12 - -C/Xii - 2^21, 



m 



= -^X21 - -E'X20, 



/ El 

— 019 020 ^10) 

P2 P2 



Y2 = 20s + {(5^ - 1)02 - 04, 
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F3 = (/5'-l)^3 + ^4 + /3^5, 



y4 = 70 - 3E^3 + Ee^ + 29j - 3105 + (1 - P^){2Ee2 - Be) + —^9, 



2p2 



021 + ~^^14 



211 + yoO'. 



5; 



Zi = (/32-l)Ci2 + 2Cii-C: 



105 



Z2 = (/52-l)Cn + C: 



10, 



Z3 — ~^Ci2 — -E'en — 2*^21! 



2(l-l/o)-y 



N2 = Ivo 



and 



TTl 1 

70 = --^^2 + E9s + -^7. 



„„ is deflned as ,„ ^ 1 + All of ,he <,ua„,i,ies C„, X™„ excepC as funcCions of the „ „ . 



come from previous work [8]. The 6*0, ... , 6*18 are found in Refs 
new functions and they are given by 

3 



9iQ, . . ., 6^22, are all 



-{E + E',){eis - e 



1 pE^, + l-p2 , r„El 

2Pp2 ' ^-^''-^ "'"^ 2^'^'' + m ' 



pi 



[3(E°)2 - /2 + 3i5(i^ + 2El)] {03 -9^- pe,) + 



P2 



{Elf 3E 3E{E + El) 1 

+ 2^2^' 2pi + 2;^^^' 



29 



- — [5(/ + ^El) + 3/?2(p22/o - /)] (^3 -e,- m - m^EEl^e, - 9s) 

+ Qp2l^yo03 + 30lE\l + PEy.o + 30El%o - ^^22 



- 6p2 



fE{p^ - 5) - 



2lpl + 2Pp2lyo{E + E'l 
h+h- 



The functions 6'i9, . . . , 6*22 are 



^12- 



with 



rp± _ 



P2 

^20 = (^20 + ^20)' 

P2 

621 = —{T21 +T21), 

P2 

622 — (^2"2+^22), 

P2 



-^20 



Xq In 



1 ± 



+ (a±)^ln 



±^0 ± a"* 



l±a± 



In 



1 ±a± 



1 ± 



- i(l T xl){l Ta^) + \{1- xl) [I - [a^f] - (1 ^ x,) [l T {a^f] 



-21 



[P2-1- Elx^ TP2yoa^ [a^li - 2] - [2a±a;o - xl + 1 + {a^J^] , 



-'-22 
21 



+ 
+ 
+ 
+ 

+ 
+ 



6± 



[{E^fxQ + p2f] - 6p2Evo{P2 + + l2mP2la^ ±E + p2 
1^^^^^=^^^! + ^ [T3(a± T l)mP2l{f3voP2 + 2E + 2E°)] A 



1 



3?7oP2 



T(p2r/o)'/5' - + I-P2T - Sp2Vo{E^J + I - P2){P2l3yo + El + ii;) h 



[E{l-p2)+P2{ElT'2l±2p2)]-3 



(«^Tl).2 



6± 



pll[±7]l + 2a^{r)o + y^)] 



It 



1 



[t/3(P277o)' - Pl{a^ T lf[^{El(3 + I-P2)- 2p2(3{a^ T 1)]] I 



^{ElP + I - p2){p2Pyo + El + E) 



{Elxtf 



6± 



(3 - 13x0 + 2/?a^) - QP2IEI 



/ ^^l/o -^0 



6± 



<^2 , 
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where = yo± a^, 6^ = 1 + /3a^, and Xq = xo + a^. The functions a^, xq, h, If, Ji, Jf, 



± T± 



and are found in Ref. Il8| 



APPENDIX C: COLLECTION OF THE pi FUNCTIONS 



Here we give the analytical results for the pi functions after performing the integrals displayed 
in Eqs. 



Pi = E'Mvo - 1) 



1-/3' 
2P 



do - Pm 



+ (E + E°)(Cio-C: 



,11 



+ P^(y,-l)9o + l'p2Y,-Plpll 



El 



Pn 



P2I 
2Mi 
1 

+ 2 



p^El^ , r^O u.. .^^nr^Oo on o„ 1 



2 + f [(2/o-l)(2i?°/3-3/)-2p2] 



r,n XlO + -7T) ^20 

2(3p2 A(3p2 



yo + 



El{El + E) 



P21 



+ 



El-E 



o "ClO + 77 

2p2 2 



Ep2 



0\2 



4p2/ 4/? 

1^ , Ei + E ^ 0p2r 



P// = 1^ [2(p2 + lyoWo - PP2012) + /(I - y'o) + 2/'r5 

E -\- E^ 
H [/3C21 + X20 - X21] 



P2i 



and 



P2I 



2Mi 



2^ 



[2^0(1 - P^)xn - 2E°xio - (1 - f3')x2i + X20] 



'^^-Cio + + ^Cii + /^P2// - ^^o' 



P2 
P2I 



P2 



E?. 



P2I 

2Mi 



4M;^° [^(^0 - 1) + 2p2] - ^X2o, 



E{p2 - /z/o)^4 + ^r/o + 1{P2 + /l/o)^5 - ^ 

P2 ^ 



+ —Cio + —AEi+mxii- 

P2 P2< 2/5p2 
f / 



E 

-X21 H 7X20 

P2I' 





E 




03- 




2 




P2 




EEl 




P2I 


Xw 



2-/5^ + 9^ 



,11 



- --(^21 - 2/^20) - 7^(i5° - E)i0i, - 210,0) + V IKvo - 1) + 2p2] 

/P2 ^P2 ^ 

All the algebraic expressions that appear in these pi and p- are defined in Appendix |Bl 
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APPENDIX D: NEW AND PREVIOUS ANALYTICAL INTEGRALS 



In this Appendix, we give a brief discussion that allows the identification of previous and new 
analytical integrals over the photon variables that emerge in the present calculation. Such integrals 
can all be put into the general form 

with a = p2, i, k and x = 1 ■ k. F and D were defined after Eq. (fT9l) . The powers of x, F, D, 
and (1 — fix) are denoted by {rpmn}. Since each one of these quantities is a rotational scalar, one 
may change the orientation of the space axes. For Si • a we use the rule jsl 

si • a ^ (si ■ p2)(a- P2), 

where a = p2, 1, or k. This considerably simplifies the calculation of such integrals. One 
can classify these integrals into three groups. In the first one we can directly identify integrals 
previously performed. This occurs for a = p2 and {rpmn} = {0012,0011,0010}, then one 
identifies 62, 63 and 64 of Ref. [8]. For a = 1 and {rpmn} = {0010,0011,0012,0121} one 
identifies the functions Cpg- For a = k and {rpmn} = {0010, 0120, 0011, 0012, 0121, 0122, 0231} 
one identifies the functions Xpq- 

In the second group are new integrals that can be expressed as combinations of previous 
results in terms of ?7o = I + yo and 6q, 62, . . . ,6is also of Ref. Jsl. This occurs for a = p2 
and {rpmn} = {0000, 0001, 00(-l)l, 0111, 0121, 0122, 0231, 1010}, for a = land {rpmn} = 
{0000,0001,0002,0120,0101}, and for a = k and {rpmn} = {0001,1010,0111,0112,0230}. 
Omitting details, such combinations are accommodated into the A/ of Appendix iBl 

The third group contains, after applying the above rule, only four new integrals with a = p2 and 
the powers are {rpmn} = {0211, 2120, 1110, 0331}. The first three of them are straightforward. 
Explicitly they are 



where 




dxx^ J dyF J dipk^ = 27r(^^2i - 2/^20), 
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and 



dxx 



-1 



yo 



-1 



2lT 



dyF I difk- = 27r(2p20 (1/0^5 - Y^] 



020 and 6^21 are new and they are listed below. 9^ is found in [8] and Yi is found in Appendix 151 
The computation of the fourth integral, 



J 



1 



yo 



2n 



although long and tedious, can be performed by using standard techniques. The final result can be 
organized as 



J 



(27r)(12P) /5 



1 



•p2?/0^3--^(/ + /?^°)^10-5E^20 



/3 

(3^ (3 



12/3 



+ P2 



P Pb+b- 
The functions 6*19, ... , 6*22 in the four new integrals are 



6*12- 



&19 
&20 
021 



xC,4{x)dx, 



x^^i{x)dx, 



-1 
1 



x'^^2ix)dx, 



and 



22 



-dx. 



1-1 1 - 

the other ^i{x), ^2{x), and ^4(0;) are used in the Ref. [6], the function ^q{x) is new and it reads 
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2 

P2V0 



1 



(X + a 



+^2 



X + a y 

3r/o [p2yo + x{El + Ix)] [El + (/ - ^2)^] 



+ 3/r/o [P2r]0 + 2x(E° + Ix)] 
1 



a 



1 



1 



+ 



0N3 
u J 



P2 



Qlpl 



X - Xo\^ - [l + xdf] 



+ 



{x + a y 

1 



+ 



X + a X + 
1 



'x + a^Y 



6+(x + a+) b^{x + a~) 



{x + a~y {x + a+y 
(1 - /3x)^4{x), 



where d = {I —p2)/El. 
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